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ABSTRACT 

A new interpretation for the second peak of T Coronae Borealis (T CrB) outbursts is proposed based 
on a thermonuclear runaway (TNR) model. The system consists of a very massive white dwarf (WD) 
with a tilting accretion disk and a lobe- filling red-giant. The first peak of the visual light curve of 
T CrB outbursts is well reproduced by the TNR model on a WD close to the Chandrasekhar mass 
(Mwd <; 1-35 Mo), while the second peak is reproduced by the combination of the irradiated M-giant 
and the irradiated tilting disk. The derived fitting parameters are the WD mass Af-yvD ~ 1-35 Ms, the 
M-giant companion mass Mrq ~ 0.7M Q (0.6 — X.OMq is acceptable), the inclination angle of the orbit 
i ~ 70°, and the tilting angle of the disk i proc ~ 35°. These parameters are consistent with the recently 
derived binary parameters of T CrB. 

Subject headings: accretion disks — binaries: close — binaries: symbiotic — novae — stars: individual 
(T CrB) 
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1. INTRODUCTION 

T Coronae Borealis (T CrB) is one of the well observed 
recurrent novae and characterized by a secondary, fainter 
maximum occurring ~ 100 days after the primary peak. 
Historically, T CrB bursted twice, in 1866 a nd 1946, wit h 
the light curves very similar each other (e.g. 
Large ellipsoidal variations in the optical ligE 
ing quiescent phase suggest that a n M3— 4 red-giant com- 
ponent fills its Roch e lobe (e.g., Leibowitz et al. 1997 ; 
Shahbaz et al. 1997| ; |Belczyriski fc Mikolajewska 1998] ). 
There have been debates on the nature of the hot compo- 
nent of this binary system. Webbink (1976) and Webbink 
et al. (1987) proposed an outburst mechanism of T CrB 
based on their main-sequence accretor model. This accre- 
tion event mode l was investigated further b y 3D numeri- 
cal si mulations ( pannizzo fc Kenyon 199*2 ; Ruffert et al 
1993]) 



However, Selvelli, Cassatella, & Gilmozzi (1992) have 
been opposed to the main-sequence accretor model from 
their analysis of IUE data in its quiescent state, which 
indicates the existence of a mass-accreting white dwarf 
(WD). The detection of X-rays and the presence of flick- 
ering in the optical light curves are also naturally ex- 
plained in terms of accretion onto a WD. Their esti- 
mated mass accretion rate in quiescent state is very high 
(M acc ~ 2.5 x 1O _8 M yr -1 ) and is exactly required by the 
thermonuclear runaway (TNR) theory to produce a TNR 
event every 80 yr on a massive (j> l.3M@) WD. Thus, the 
1866 and 1946 outbursts can be interpreted in terms of a 
TNR event on a very massive WD. 

Rapid decline rates of the light curves indicate a very 
massive WD clo se to the C handrasekhar limit, J1%d ~ 
1.37 - 1.38M ( |Kato 1995| , 1999). Assuming solar com- 
position of the WD envelope, Kato calculated nova light 



curves for WD masses of 1.2, 1.3, 1.35 and 1.377 M Q and 
found that the light curve of the 1.377M Q model is in bet- 
ter agreement with the observational light curve of T CrB 
than the other lower mass models. 

Recently, other observational supports for a massive 
WD in T CrB have been reported. Belczyhski and Mikola- 
jewska (1998) derived a permitted range of binary param- 
eters, AfwD = 1.2 ± 0.2Af Q , from amplitude of the ellip- 
soidal variability and constraints from the orbital solution 
of M-giants. In Shahbaz et al. (1997), a massive WD of 
Mwd = 1-3 — 2.5M Q is suggested from the infrared light 
curve fitting. Combining these two permitted ranges of the 
WD mass in T CrB, we may conclude that a mass of the 
WD is between Mwd = L3 — 1.4M©, which is very consis- 
tent with the light curve analysis Mwd ~ 1-37 — 1.38M 
by Kato (1999). 

The secondary maximum in outbursts is not generally 
observed in classical novae or in other recurrent novae. 
Selvelli et al. (1992) suggested a possibility of irradiation 
by a stationary shell around the system, although the pres- 
ence of the shell is just a speculation. In this Letter, we 
propose another possibility of the second peak: irradiation 
by a tilting accretion disk around a massive WD. 

The main results of our analysis are: 1) the first peak is 
naturally reproduced by the fast developing photosphere 
of the WD envelope based on the TNR model incorporated 
with a very massive WD (Mwd ~ 1.35M©). 2) the sec- 
ond peak is not fully reproduced by an irradiated M-giant 
model as simply estimated by Webbink et al. (1987); 3) 
the second peak can be well reproduced if we introduce 
an irradiated tilting accretion disk around the WD in ad- 
dition to the irradiated M-giant companion. Such tilting 
instabilities of an accretion disk have been suggested by 
Pringle (1996) for central stars as luminous as the Edding- 
ton limit or more (radiation- induced instability). Because 
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the maximum luminos ity of T CrB outbu rsts exceeds the 



Eddington limit (e.g., Selvelli ct al. 1992), the radiation- 



induced instability may work well in T CrB outbursts. 

2. THEORETICAL LIGHT CURVES 

Our model is graphically shown in Figure H. The visual 
light is contributed to by three components of the system: 
1) the WD photosphere, 2) the M-giant photosphere, and 
3) the accretion disk surface. 



2.1. Decay Phase of Novae 



In the TNR model, WD envelopes expand greatly as 
large as ~ 100 i?© or more and then the photospheric 
radius gradually shrinks to the original size of the white 
dwarfs (e.g., ~ 0.004 R Q for M WD = 1.35 M Q ). The 
optical luminosity reaches its maximum at the maximum 
expansion of the photosphere and then gradually darkens 
to the level in quiescent phase. Since the WD enve lopes 
reach a steady-st ate in the decay phase of novae (e.g. , Katcj 
& Hachisu 1994), we are able to treat the development of 
the envelope by a unique sequence of steady-state solu- 
tions with different envelope masses (AM) as shown by 
Kato & Hachisu (1994). 

We have calculated such sequences for WDs with various 
masses of M WD = 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.35 
and 1.377M© and obtained the optical light curves for the 
decay phase of TNR events. Here, we choose 1.377M© as 
a limiting mass ju st below the mass at the SN la explosion 
in W7 (1.378M©, |Nomoto ct al. 1984|). We have u sed the 
updated OPAL opacity ( |[glesias fc Rogers" 1996| ), which 
has a strong peak near logT ~ 5.2 abou t 20 — 30% larger 



than that of the original OPAL opacity ( Rogers & Iglesias 
1992[ ) which was used in Kato & Hachisu (1994). The nu- 
merical method and various assumptions are the same as 
those in Kato & Hachisu (1994). It should be noted here 
that optically thick winds blow when the WD envelope ex- 
pands and the photospheric temperature decreases below 
logT ph ~5.5. 

Each wind solution is a unique function of the enve- 
lope mass AM if the WD mass is given. The envelope 
mass is decreasing due to the wind mass loss at a rate 
of M w i n( j(AM) and hydrogen shell burning at a rate of 
M nuc (AM), i.e., 



— AM = M acc 
at 



Mv ind - M n 



(1) 



where M acc is the mass accretion rate to the WD. Inte- 
grating equation (|j), we follow the development of the 
envelope mass AM and obtain physical quantities such 
as the photospheric temperature T p h, photospheric radius 
R p h, photospheric velocity v p h, wind mass loss rate M w ; n d, 
and nuclear burning rate M nuc . When the envelope mass 
decreases to below the critical mass, the wind stops and 
after that the envelope mass is decreased only by nuclear 
burning. 

2.2. White Dwarf Photosphere 

We have assumed a black-body photosphere of the white 
dwarf envelope. After the optical peak, the photosphere 
shrinks with the envelope mass being blown off in the 
wind, and the photospheric temperature increases with 



the visual light decrease because the main emitting region 
moves blueward. Based on our wind solutions, we have 
obtained visual magnitude of the WD photosphere with a 
window function given by Allen (1973). The photospheric 
surface is divided into 16 pieces in the latitudinal angle 
(A9 = 7r/16) and into 32 pieces in the longitudinal angle 
(A</> = 27r/32) as shown in Figure [l]. Then, the contribu- 
tion of each piece is summed up by considering the incli- 
nation angle to the viewer. A linear limb-darkening law 



(the coefficient is x = 0.95, see Belczyhski & Mikolajewska 



1998; Wilson 1990) is incorporated into the calculation. 



2.3. Companion M-giant Photosphere 

To construct a light curve, we have also included the 
contribution of the companion star irradiated by the WD 
photosphere. The surface of the companion star is as- 
sumed to fill the inner critical Roche lobe as shown in 
Figure ^. Numerically dividing the latitudinal angle into 
32 pieces (A8 — tt/32) and the longitudinal angle into 64 
pieces (A<f> = 27r/64), we have also summed up the con- 
tribution of each area considering the inclination angle to 
the viewer by assuming the same limb-darkening law as the 
WD photosphere, but we neglect the gravity-darkening ef- 
fect of the companion star because the hemisphere to the 
WD is heated up by the irradiation. Here, we assume that 
50% of the absorbed energy is reemitted from the hemi- 
sphere of the companion with a black-body spectrum at 
a local temperature. The original (non-irradiated) photo- 
spheric temperature o f the companion star is assumed t o 
be Tph.RG = 3500 K (|Bclczyhski fc Mikolajewska 1998| ). 
If the accretion disk around the WD blocks the light from 
the WD photosphere, it makes a shadow on the surface of 
the companion star. Such an effect is also included in our 
calculation. 

The orbit of the companion star is assumed to be cir- 
cular. The light curves are calculated for five cases of 
the companion mass, i.e., Mrg = 0.6, 0.7, 0.8, 0.9, and 
l.OAf©, as derived by Belczynski & Mikolajewska (1998). 
Since we obtain similar light curves for all of these five 
masses, we show here only the results for Mrg = 0.7M©. 
In this case, the separation is a = 199. 3i?©, the effective 
radius of the inner critical Roche lobe for the WD compo- 
nent is R\ = 87.Oi?0, the effective radius of the M-giant 
companion star is i?2 = R2 — 64. 5i?©. 

2.4. Accretion Disk Surface 

We have included the luminosity coming from the ac- 
cretion disk irradiated by the WD photosphere when the 
accretion disk reappears a few to several days after the op- 
tical maximum. The surface of the accretion disk absorbs 
photons and reemits in the same way as the companion 
does. Here, we assume that the accretion disk surface 
emits photons as a black-body at a local temperature of 
the heated surface. We assume further that 25% of the ab- 
sorbed energy is emitted from the surface, while the other 
is carried into interior of the accretion disk and eventually 
brought into the WD. The temperature of the outer edge 
is assumed to be Td; s k = 2000 K, which is never irradiated 
by the WD photosphere. 

An axisymmetric accretion disk with a thickness given 

by 

h = /3R disk ( -^-) , (2) 
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is assumed, where h is the height of the surface from the 
equatorial plane, w the distance on the equatorial plane 
from the center of the WD, i?disk the outer edge of the 
accretion disk, and (3 is a numerical factor showing the 
degree of thickness. Here, we assume (3 = 0.01 during the 
strong wind phase because the flaring-up edge of the accre- 
tion disk is blown in t he wind, but it incre ases to (3 = 0.15 
after the wind stops ( [Bchandl et al. 1997 ). The surface of 
the accretion disk is divided into 16 pieces logarithmically 
in the radial direction and into 32_pieces evenly in the az- 
imuthal angle as shown in Figure [II The outer edge of the 
accretion disk is also divided into 32 pieces by rectangles. 

The luminosity of the accretion disk depends strongly 
on both the thickness of f3 and the radius of i?disk- The 
size of the accretion disk is also assumed to be given by 



-Rdisk = aR 1 , 



(3) 



where Rl is the effective radius of the inner critical Roche 
lobe given by Eggleton's (1983) formula. The viscous heat- 
ing is neglected because it is rather smaller than that of 
the irradiation effects. 

We assume that the disk is tilting due to Pringlc's (1996) 
mechanism. The tilting disk is introduced by inclining the 
above disk by degree of z prec = 35° with a precessing an- 
gular velocity of 

^prcc T^orb- (4) 

The initial phase of precession is assumed to be 0o = 
— 170° at the epoch of spectroscopic conjunction with the 
M-giant in front, i .e., JD 243 1931.05 + 227.67-B at E = 
flLincs ct al. 198S| ). 



3. RESULTS 

To fit the first peak of the outburst light curve, we have 
calculated four cases of T^-magnitude light curves with four 
different WD masses, i.e., M WD = 1.377 M , 1.35 M Q , 
1.3 M Q , and 1.2 M Q and found that both the 1.377M Q and 
1.35M Q light curves are in much better agreement with the 
observational one than the other less massive ones as al- 
ready shown by Kato (1995, 1999). Therefore, we adopt 
the Mwd = 1-35 M Q in this Letter. 

The optically thick wind stops about 100 days after the 
luminosity peak, around 1960 (JD 2430000+) as shown in 
Figure 0. Here, we assume that the mass accretion rate 
(Mace) to the WD is increased to a few times 10~ 7 Mq 
yr _1 , because the disk is heated up by the WD photo- 
sphere (also once engulfed by the WD photosphere). The 
hydrogen shell-burning vanishes 120 days after the peak 
when M acc < 1 x 10~ 7 M Q yr _1 , while it continues for 
rather long period when M acc > 4x 10~ 7 A/q yr _1 as shown 
in Table |f 

The upper panel in Figure |^ shows the theoretical V- 
magnitude light curve (solid line) of the WD photosphere 
with M acc = 4 x lO~ 7 ill0 yr _1 and the companion without 
irradiation together with the ob servational points (open 
circles; taken from Pettit 1946a, b, c, d). Two arrows 



indicate epochs o f the spectroscopic conjunction with the 
M-giant in front (Lines et al. 1988). The ellipsoidal light 
variation is clearly shown in the later phase. The second 
panel depicts the T^-magnitude including the effects of the 
companion irradiated by the WD. As already suggested 



by Webbink et al. (1987), this cannot well reproduce the 
second peak of the light curves. The third panel shows the 
theoretical light curve further including the tilting accre- 
tion disk irradiated by the WD photosphere. It fully repro- 
duces the observational light curve if we choose the fitting 
parameters of a — 0.7, (3 = 0.15, 7 = 9/8, i pl - cc = 35°, and 
O = -170°. 

4. DISCUSSION 

The radiation induced instability of the accretion disk 
sets in if the condition 
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(5) 



is satisfied flSouthwcll ct al. 1997j ). Selvelli et al. (1992) 
estimated the accretion rate of T CrB as M acc ~ 0.25 x 
10~ 7 Mq yr _1 , which meets condition (0). Therefore, we 
can expect the radiation induced instability in T CrB sys- 
tem. The growth timescale of warping is estimated by 
Pringle (1996) and Livio & Pringle (1996) as the same 
timescale of the precessing period, i.e., 
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A/disk 



1 x 10" 6 M 
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2 x 10 38 erg s" 1 



day , (6) 



where we assume that the a-parameter of the standard 
accretion disk is a ~ 0.1. Thus, the growth timescale is 
short enough to excite warping of the accretion disk. It 
should be noted here that the heated accretion disk sur- 
face, just after the maximum expansion of the WD photo- 
sphere, could drive a disk wind because the disk has once 
been engulfed by the WD photosphere. Such a wind could 
exert an even larger back pressure on the disk than radia- 
tion and thus lead to tilting (e.g., |Schandl fc Meyer 1994 ; 
Schandl 1996| ). Its growth timescale is much shorter than 
the timescale given by equation @. 

The large inclination angle of the accretion disk such 
as i proc ~ 30 — 35° is required to reproduce the obser- 
vational light curve mainly because we need a reflection 
area as large as the companion, which is viewed from ours 
as shown in Figure ||. About 10% faster precessing an- 
gular velocity is also required from the phase relation be- 
tween the rising shoulder of the second peak near 1970 (JD 
2430000+) , a small dip near 2100 JD, and then a small 
bump near 2130 JD as shown in Figure |. These dips are 
caused by a large shadow on the companion cast by the 
accretion disk. This precession velocity is consistent with 
the fact that the radiation-induced precession is prograde, 
while the tidally induced precession is retrograde. 
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Sports. 
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Table 1 

Periods of wind and nuclear burning phases after the peak. 
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wind phase 
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93 
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Fig. 1. — Model configuration near the second peak of the recurrent nova T CrB. The cool component (right in the figure) is a red-giant 
filling up its inner critical Roche lobe. The hemisphere is heated up by the hot component (1.35 Mq white dwarf, left in the figure). We 
assume a tilting accretion disk around the hot component, which is precessing about 10% faster than the orbital rotation. The surface of the 
accretion disk is also heated up by the hot component. The photosphcric radius of the hot component near the second peak is as small as 
~ 0.05 Rq, about ~ 0.001 times the size of the cool component, but it is exaggerated in this figure to easily see it. 
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Fig. 2. — Model light curves are plotted against time (JD 2430000+) together with the observational points (Pettit 1946a, 1946b, 1946c, 
1946d). Top panel: the visual magnitude of the sum of the white dwarf (WD) photosphere and the red giant (RG) photosphere without 
irradiation. Large arrows indicate epochs at the spectroscopic conjunction with the M-giant in front. Middle panel: the visual magnitude of 
the WD photosphere and the RG photosphere with irradiation. Bottom panel: the visual magnitude of the WD, RG with irradiation, and 
the accretion disk hcatcd-up by the hot component. 



